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Abstract. We present a compilation of PRONAOS-based results concerning the temperature dependence of the 
dust submillimeter spectral index, including data from Galactic cirrus, star-forming regions, dust associated to a 
young stellar object, and a spiral galaxy. We observe large variations of the spectral index (from 0.8 to 2.4) in a 
wide range of temperatures (11 to 80 K). These spectral index variations follow a hyperbolic-shaped function of 
the temperature, high spectral indices (1.6-2.4) being observed in cold regions (11-20 K) while low indices (0.8-1.6) 
are observed in warm regions (35-80 K). Three distinct effects may play a role in this temperature dependence: 
one is that the grain sizes change in dense environments, another is that the chemical composition of the grains 
is not the same in different environments, a third one is that there is an intrinsic dependence of the dust spectral 
index on the temperature due to quantum processes. This last effect is backed up by laboratory measurements 
and could be the dominant one. 
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1. Introduction 



Characterizing dust emissivity properties is one of the ma- 
jor challenges of nowadays submillimeter astronomy. In 
this spectral range, big grains at thermal equilibrium (see 
for instance De sert et al. 1 9901 dominate the dust emis- 
sion. Their emission is characterized by a temperature and 
a spectral dependence of the emissivity which is usually 
simply modelled by a spectral index. 

The temperature, density and opacity of a molecular 
cloud are key parameters which control the structure and 
evolution of the clumps, and therefore, star formation. The 
spectral index (/3) of a given dust grain population is di- 
rectly linked to the internal physical mechanisms and the 
chemical nature of the grains. 

It is generally admitted from Kramcrs-Konig re- 
lations that 1 is a lower limit for the spectral index 
IjEmerson 1988|l . A classical value is 2, following the inves- 
tigations of Gezari et al. (1973) and the calculations of, for 
instance, Draine & Lee (1984). This value is particularly 
invoked for isotropic crystalline grains, and it is thought 



infrared: ISM: continuum to be an upper limit IjWooten 1972L lEmerson 1988|l . 

For amorphous silicate or graphitic grains, a value of 
P — 2 is also favoured for different physical reasons 
URowan-Robinson 19861 ITielens fc Allamandola 19871 
IMennella et al. 1995|l . However, it is not the case for 
amorphous carbon, which is thought to have a spectral 
index equal to 1 IjKoike et al. 1980)l . as well as aggre- 
gates of silicates and graphite in a porous structure 
CMathis fc Whiffen 198911 . Silicate grains embedded in ice 
mantles are thought to have spectral indices between 1.5 
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and 2 |Miyake fc Nakaga wa 1993[ IPreibisch et al. 19931 
IKriigel fc Siebenmorgcn 1994[ IPollack et al. 1994|l . In 
general, amorphous rather than crystalline structure and 
increase of the grain size ( [Mannings fc Emerson 1994| ) 
are the main reasons to predict submillimeter spectral 
indices lower than 2. Spectral indices above 2 may 
exist, according to several laboratory measurements on 
grain analogs IjKoike et al. 19951 IMennella et al. 1995L 
jAgladze et al. 1996| IMennella et al. 1998|l . Also, Agladze 
et al. (1996) and Mennella et al. (1998) showed an 
anticorrelation existing for some types of grains between 
the temperature and the emissivity spectral index. 

Observations of the diffuse interstellar medium 
at large scales favour /3 around 2 (see for exam- 
ple IBoulanger et al. 1996| IDunne fc Eales 2001|l . In the 
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case of molecular clouds, spectral indices are usually 
found to be between 1.5 IjWalker et al. 1990)l and 2 
(e.g. [Wright et al. 1992] in Orion). However, low val- 
ues (0.2-1.4) of the spectral index have been observed 
in circumstellar environments HWeintraub et al. 19891 
IBeckwith et al. 19901 [knapp et al. 1993| ). This can be at- 
tributed to grain growth in dense stellar envelopes. Low 
indices have also been observed in molecular cloud cores 
fe.g-. IBlake et al. 1996|l . Spectral indices larger than 2 have 
also been observed in the millimeter range l|Schwartz 19821 
IWilson et al. 19951 IKuan et al. 1996|l . in particular on 
the whole sky at large scales by the WMAP satellite 
HBennett et al. 2003|l . 

In this context, submillimeter multi-band imaging of 
the interstellar medium provides a lot of useful informa- 
tion about dust properties and interstellar medium struc- 
ture, especially if the dust emission parameters, namely 
both the temperature and the spectral index, can be prop- 
erly derived on top of submillimeter intensities. 

2. Observations 

PRONAOS (PROgramme NAtional d'Observations 
Submillimetriques) is a French balloon-borne submillime- 
ter experiment IjSerra et al. 20011 IRistorcelli et al. 19981 
ILamarre et al. 1994|l . Four bolometers cooled at 0.3 
K measure the submillimeter flux with sensitivity to 
low brightness gradients of about 4 MJy/sr in band 
1 (200 urn) and 0.8 MJy/sr in band 4 (580 ^im). The 
effective wavelengths are 200, 260, 360 and 580 /xm, 
and the angular resolutions are 2' in bands 1 and 2, 
2.5' in band 3 and 3.5' in band 4. The data that we 
analyze here were obtained during the second flight 
of PRONAOS in September 1996, at Fort Sumner, 
New Mexico. The data processing method, including 
deconvolution from chopped data, is described in Dupac 
et al. (2001), and the calibration procedure is detailed 
in Pajot et al. (2003). This experiment has observed 
various phases of the interstellar medium, from diffuse 
clouds in Polaris ([Bernard et al. 1999|l and Taurus 
l |Stepnik ct al. 2003) to massive star-forming regions in 
Orion (RistorceUi et al. 1998. and [Dupac et al. 20011 ), 
Messier 17 ( [Dupac et al. 2002[ ), Cygnus B, and the 
dusty envelope surrounding the young massive star 
GH20 092.67-1-03.07 in NCS IRernard et al. 200311 . The 
p Ophiuchi low-mass star-forming region has also been 
observed, as well as the edge-on spiral galaxy NGC 891 
l [Dupac et al. 2003| ). In this letter, we analyze constraints 
that PRONAOS data establish on the dust spectral index 
and its relation to the temperature. Ristorcelli et al. 
(1998) and Bernard et al. (1999) gave tentative evidences 
of a spectral index larger than 2, respectively in Orion 
and the Polaris flare. Dupac et al. (2001) and Dupac 
et al. (2002) showed large variations of the spectral 
index, including values above 2, and the presence of an 
anticorrelation between the temperature and the spectral 
index in massive star-forming regions. In this letter. 



we extend the analysis to other regions observed with 
PRONAOS. 

3. Analysis 

3.1. Bases of the analysis 

For Orion ( [Dupac et al. 2001D, M17 ( [Dupac et al. 2002[ ) 

and Cvgnus ( Boud et et al. 2003jl , we analyze the intensity 
maps using the following procedure. With the PRONAOS 
and IRAS intensity maps, we produce temperature (T) 
and spectral index {(3) maps following the processing de- 
scribed in Dupac et al. 2002 (see also [Dupac et al. 200 1| ). 
We fit a modified black body law to the spectra: Ii, = 
£0 By{\,T) {X/Xo)~^, where Ii, is the spectral intensity 
(MJy/sr), eo is the emissivity at Aq of the observed dust 
column density, B^, is the Planck function, T is the tem- 
perature and P is the spectral index. In most of the ar- 
eas, we use either only PRONAOS data or PRONAOS -I- 
IRAS 100 p,m data. The reasons why we do not use IRAS 
data everywhere are that first, some very intense regions 
are saturated in the 100 fim IRAS band, second, many 
faint regions which are clearly detected with PRONAOS 
are very noisy with IRAS, or a large-scale warm compo- 
nent dominates in the IRAS bands. In the non-saturated 
warmest areas (> 70 K), we use IRAS 60 ^m data too 
(see [Dupac et al. 2002[ for more details). 

In this work, we restrain the analysis to all fully inde- 
pendent (3.5' side) pixels for which both relative errors on 
the temperature and the spectral index are less than 20%. 
This procedure is motivated notably by the fact that there 
exists a degeneracy between the temperature and the spec- 
tral index. This effect is mainly due to the relatively large 
error bars on the spectral index in cold regions and on 
the temperature in warm regions (Rayleigh- Jeans limit). 
Dupac et al. (2001) and Dupac et al. (2002) have shown by 
fitting simulated data that this artificial anticorrelation ef- 
fect was small compared to the effect observed in the data. 
For this work, we have also done other simulations, in a 
slightly different way from the previous ones: we consider 
temperature values in the same temperature range as the 
observations, but with a constant spectral index (1.5), we 
produce simulated data with a realistic noise level, and 
we fit them with the usual procedure. With PRONAOS 
simulated data alone, we find that the relative residual er- 
ror on the spectral index is 22% rms. With PRONAOS + 
IRAS 100 /im data, the same error is only 7% rms on the 
spectral index and 9% on the temperature. The induced 
inverse correlation effect is very limited. Thus, this noise 
effect is clearly not capable to explain the T-f] distribution 
found in the data. Here we still reduce it by restraining 
the analysis to pixels for which errors on both the tem- 
perature and the spectral index are lower than 20%. Also, 
we disregard pixels for which at least one of the bands has 
a spectral intensity below 10 MJy/sr. Moreover, all spec- 
tra were visually checked in order to exclude pixels which 
could be contaminated by noise features. This provides 
us with a set of 31 pixels for M17, 19 for Orion and 60 
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Wovelenqth (/^m) 

Fig. 1. Examples of spectra. From top to bottom, 
the crosses represent data from the Orion Molecular 
Cloud - 1, M17 North, and Cloud 2 in Orion (see 
|Dupac et al. 2001] for details). Triangles represent data 
from the south-western region of the NGC 891 galaxy (see 
|Dupac et al. 2003| ). The error bars are plotted within the 
3 a confidence level. The full lines represent the results 
of the fit by a single modified black body. The parame- 
ters of each fit (temperature and spectral index) are also 
indicated. 



for Cygnus. For the cirrus clouds in Polaris and Taurus, 
we have, respectively, two and one data points avail- 
able HBernard et al. 1999l and |Stepnik et al. 2003| ). For the 
NCS data point, we fit the PRONAOS fluxes as given in 
Bernard et al. (2003). We also present five data points 
from p Ophiuchi IjRistorcelli et al. 2003|l and three from 
the NGC 891 edge-on spiral galaxy ( |Dupac et al. 2003| ). 
We present typical examples of the obtained spectra in 
Fig. ID 



3.2. Dependence on the temperature of the spectral 
index 

We present in Fig. 13 and 13 temperature versus spectral 
index plots of the above-defined pixel set. For clarity, we 
present all points in Fig. |21 without the error bars, and a 
plot with the error bars in Fig. |2] but without the numer- 
ous points from Cygnus and M17. These data points have 
similar error bars to the Orion data points. The temper- 
ature in this data set ranges from 11 to 80 K, and the 
spectral index also exhibits large variations from 0.8 to 
2.4. One can observe an anticorrelation on these plots be- 
tween the temperature and the spectral index, in the sense 
that the cold regions have high spectral indices around 2, 
and warmer regions have spectral indices below 1.5. In 
particular, no data points with T > 35 K and (3 > 1.6 can 
be found, nor points with T < 20 K and (3 < 1.5. All cold 
regions below 18 K have spectral indices above 1.8, except 
one point with a large error bar, which is well compatible 
with a value above 1.8. 




Temperature (K) 

Fig. 2. Spectral index versus temperature, for fully in- 
dependent pixels in Orion (black), p Ophiuchi (green), 
Polaris (blue), Taurus (purple), NCS (hght blue) and NGC 
891 (red). The error bars are plotted within the 68% con- 
fidence interval. For clarity, the M17 and Cygnus points 
are not presented in this figure, but they have similar error 
bars as the Orion points. 




Fig. 3. Spectral index versus temperature, for fully inde- 
pendent pixels in Orion (red asterisks), M17 (black di- 
amonds), Cygnus (orange triangles), p Ophiuchi (green 
crosses), Polaris (blue squares), Taurus (purple square), 
NCS (hght blue cross) and NGC 891 (purple crosses). 
The full line is the result of the best hyperbolic fit: 
B = 1 

1-^ 0.4+0.008T 



This anticorrelation effect is present for all objects in 
which we observe a large range of temperatures, namely 
Orion, M17, Cygnus and p Ophiuchi. It is also remarkable 
that the few points from other regions are well compati- 
ble with this general anticorrelation trend. Following the 
type of the observed object, the temperature distribution 
of the points is different, but the points seem to follow 
this same general trend of an inverse dependence of the 
spectral index on the temperature. Hereafter, we therefore 
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consider that this effect can be simply investigated with 
the different objects altogether. 

Fitting the data under the assumption of a constant 
spectral index leads to / degrees of freedom (d.o.f. 
hence on) values of 1835/121 for /3=2 and 377/121 for 
P=1.5. These high values of the x^/d.o.f. exclude the pos- 
sibihty of matching the data with an invariant spectral 
index. The hnear regression fit of P{T), considering the 
temperature as an independent variable, gives a slope of 
-0.014 K-\ and /3 = 2.01 at T = 0. Then the linear cor- 
relation coefficient is —0.79, and the x^/d.o.f. is 137/120 
(=1-14). 

Since the points seem quite well distributed along a 
hyperbolic-shape line, we perform a hyperbohc fit of the 
data using: /? = j^^, where 6 and oj are free parame- 
ters. Considering that the temperature is the independent 
variable, the best fit gives (5=0.40 ± 0.02 and u;=0.0079 
± 0.0005 K-i, with xVd-o.f. = 120/120. The tempera- 
ture dependence of the emissivity spectral index is thus 
very well fitted with an hyperbolic approximating func- 
tion. Note that since the number of free parameters is the 
same as in the linear regression, the lower x^ is significant, 
so that we can consider the hyperbolic fit to be better. 

If now we consider that the spectral index is the inde- 
pendent variable, fitting the T(/3) function by the linear 
regression gives a large x^/d.o.f. value of 1099/120. In this 
case, the hyperbolic fit gives x^/d.o.f. = 1123/120. These 
bad fits are mainly due to the large dispersion of the spec- 
tral indices that we observe at low temperatures. 

4. Discussion 

In order to interpret this inverse temperature dependence 
of the spectral index, several hypotheses can be investi- 
gated. Mixtures of dust components with different temper- 
atures, but having the same spectral index, can be put for- 
ward to explain apparent low indices. However, this effect 
was shown to be small compared to the amplitude of the 
observed spectral index variations ( |Dupac et al. 200'2| ). 

As stated in the introduction, low indices have been ob- 
served in active environments such as circumstellar disks 
and warm molecular cloud cores. PRONAOS observations 
of Orion, M17, Cygnus B and p Ophiuchi show such low 
spectral indices around 1. This could be due to the increase 
of the grain size in dense environments, as predicted by 
Miyake & Nakagawa (1993), whose models with different 
grain size distributions show that the maximum grain ra- 
dius adopted is crucial for the long-wavelength emissivity 
slope. This is also the reason invoked by Goldsmith et al. 
(1997) to explain their results, which showed an inverse 
dependence of the spectral index on the optical depth. 
However, we do not find such correlation in our data. 
Moreover, the range of variations that they find for the 
spectral index (from 1.6 to 3.2) is significantly different 
from ours (from 1 to 2.4), though with a similar variation 
amplitude. Also, these grain growth effects are invoked for 
very peculiar dense environments around stars rather than 
for large scale observations such as the PRONAOS ones. 



Therefore, it is difficult to consider that size distribution 
effects are dominant in the T-P anticorrelation that we 
observe, though they may play a role. 

Another slightly different explanation could be that 
one finds different chemical composition or physical 
state of the grains following the physical conditions 
of the medium, in particular following the temper- 
ature. Cold environments could harbor fluffy silicate 
grains including ice compounds, having spectral indices 
of 2 (e.g. IKriigel fc Siebenmorgen 1994|), or simply sili- 
cate or graphitic particles (e.g. Ibraine &: Lee 1984|l . or 
even olivine core or fused quartz grains with ice man- 
tles IjAannestad 19751 /? « 3) to explain high indices 
above 2. Warm regions could harbor aggregates of sili- 
cates, porous graphite or amorphous carbon and therefore 
have a spectral index around 1 (|Mathis fc Whiffen 1989|l . 
Nevertheless, more observations are needed to constrain 
the chemical structure of the large grains in various me- 
dia and confirm this quite speculative hypothesis. 

Finally, an explanation may be searched for in a 
temperature dependence of the intrinsic optical prop- 
erties of the grains (see [Henning fc Mutschke 1997| for 
a review of the different physical processes which 
could be involved). This can be linked to laboratory 
results dAgladze et al. 1996| and IMennella et al. 1998|l 
which showed a T-/3 anticorrelation on interstellar grain 
analogs. Mennella et al. (1998) measured the absorption 
coefficient of crystalline and amorphous grain analogs be- 
tween 20 /im and 2 mm wavelength, in the temperature 
range 24-295 K. They derived an anticorrelation between 
the temperature and the spectral index, and attributed 
it to two-phonon difference processes ( [Sparks et al. 1982| ) . 
Agladze et al. (1996) measured absorption spectra of sil- 
icates between 0.7 and 2.9 mm wavelength and deduced 
an anticorrelation between the spectral index and the tem- 
perature in the temperature range 10-25 K, for two dif- 
ferent silicate precursors. They attributed the underly- 
ing absorption process to a resonant tunneling effect be- 
tween ground states of two-level systems ( [Phillips 1972| ). 
The temperature dependence of the spectral index which 
they show is in a range of temperatures and spectral 
indices close to our measurements: they measure spec- 
tral indices between 1.8 and 2.6 for temperatures rang- 
ing from 10 to 25 K, while we measure spectral in- 
dices roughly between 1.4 and 2.4 in the same tem- 
perature domain. However, the wavelength range is not 
the same, and there is no evidence that this effect is 
simply transposable from the millimeter range to the 
PRONAOS submillimeter range, since experimental re- 
sults suggest a cut-off around 500 /im for this phenomenon 
f Fitzgeral d et al. 2000| ). Furthermore, at submillimeter 
wavelengths and higher temperatures, thermally-activated 
relaxation processes and temperature-dependent absorp- 
tion associated to transitions to excited levels of two-level 
systems ( jFitzgerald et al. 2000i IHutt et al. 1989|l might 
dominate. These issues will be treated in a forthcoming 
paper (Meny et al., in prep.) This possible discrepancy 
between the submillimeter and the millimeter spectral in- 
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dices could explain the small difference observed between 
the WMAP (millimeter) spectral index of 2.2 and the 
PRONAOS one, which is around 2 at DIRBE (rather low) 
large-scale temperatures. 

5. Conclusion 

We have presented a compilation of PRONAOS-based re- 
sults concerning the temperatures and spectral indices 
measured in a variety of regions of the interstellar medium. 
The temperatures observed range from 11 to 80 K, which 
provides a large field of investigation for measuring vari- 
ations of the dust emissivity. We give strong evidence of 
large variations of the submillimeter spectral index be- 
tween 0.8 and 2.4, as well as of the existence of an in- 
verse dependence of the emissivity spectral index on the 
temperature. This dependence is well fitted by the func- 
tion j3 = Q ^j^Q QQg/p ■ Several interpretations are possible 
for this effect: one is that the grain sizes change in dense 
environments, another is that the chemical composition 
of the grains is not the same in different environments 
and that this correlates to the temperature, a third one 
is that there is an intrinsic dependence of the spectral in- 
dex on the temperature, due to quantum processes such as 
two-level tunneling effects. Additional modeling, as well as 
additional laboratory measurements and astrophysical ob- 
servations, are required in order to discriminate between 
these different interpretations. However, the most likely 
explanation could be that there is an intrinsic dependence 
of the emissivity on the temperature. A new dust semi- 
empirical model is needed to take into account these vari- 
ations of the spectral index. 
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